A toy model of fractal glioma development under RF electric field treatment 
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A toy model for glioma treatment by a radio frequency electric field is suggested. This low- 
intensity, intermediate-frequency alternating electric field is known as the tumor-treating-field 
(TTF). In the framework of this model the efficiency of this TTF is estimated, and the inter- 
play between the TTF and the migration-proliferation dichotomy of cancer cells is considered. The 
model is based on a modification of a comb model for cancer cells, where the migration-proliferation 
dichotomy becomes naturally apparent. Considering glioma cancer as a fractal dielectric composite 
of cancer cells and normal tissue cells, a new effective mechanism of glioma treatment is suggested 
in the form of a giant enhancement of the TTF. This leads to the irreversible electroporation that 
may be an effective non-invasive method of treating brain cancer. 

PACS numbers: 05.40.FbRandom walks and Levy flights, 87.50.S-Radiofrequency/microwave fields effects 
, 87.17.EeGrowth and division 



I. INTRODUCTION 

Recent progress in treating glioma using an alternat- 
ing, radio frequency electric field opens new questions 
about the migration-proliferation dichotomy of cancer 
cells extended to RF electric field. It was shown in 
in vitro and in vivo experiments that a low-intensity, 
intermediate-frequency (100-300 kHz), alternating elec- 
tric field, known as the tumor- treating field (TTF), de- 
stroys cells that are undergoing division [lJ-Q ■ As ex- 
plained, the modality of the TTF treatment is physical 
and related to dielectrophoresis (see e.g., [4|), when the 
non-uniformity of the TTF exerts a force, focusing at the 
narrow cytoplasmatic bridge between two daughter cells 
at fission. This leads to disruption of the cell prolifera- 
tion, and, eventually, the cells are destroyed without the 
quiescent cells of normal tissues being affected [l|-y| . 

One of the main features of malignant brain cancer 
is the ability of tumor cells to invade the normal tis- 
sue away from the multi-cell tumor core, and this motil- 
ity is the most critical feature of brain cancer, causing 
treatment failure J5J. Clinical investigations of glioma 
show that the proliferation rate of migratory cells is es- 
sentially lower in the invasion region than in the tumor 
core [a, |7J . This phenomenon is known as the migration- 
proliferation dichotomy, which is an inherent property of 
glioblastoma (glioma) , making it a highly aggressive and 
invasive tumor. A question that may be of great impor- 
tance for glioma treatment is whether the TTF is effective 
against invasive cells in the presence of the migration- 
proliferation dichotomy, when switching between migrat- 
ing and proliferating phenotypes takes place. The switch- 
ing process between these two phenotypes is still not well 
understood, and relevant models with different switch- 
ing mechanisms are suggested. Comprehensive reviews 
of these sophisticated approaches can be found [8h12|. 
Therefore, for the medical treatment of cancers, dielec- 
trophoresis of live cancer cells can be a new direction 
towards understanding this switching, since the polar- 
ization properties of migrating and dividing cells differ 
[lj, [l5|, and this great selectivity can be manipulated. 



One of these experimental possibilities is dielectric spec- 
troscopy, which is a powerful tool for investigating the 
dielectric properties of cells and the structural parts of 
cells, and can provide important knowledge about differ- 
ent cell structures, and resolve the different properties of 
normal and malignant cells [lj, [l5[ . 

Therefore, an essential question is how the TTF affects 
aggressive migrating cells in the outer-invasive region 
with a low-rate of proliferation, and what should be done 
to make the TTF an effective means of manipulating the 
cancer in this region. To shed light on this situation, we 
consider a simplified toy model, where a mathematical 
formulation of the migration-proliferation dichotomy can 
be performed in the framework of continuous-time ran- 
dom walks (CTRW) [MUzl- Tne simplest mathematical 
realization of the CTRW mechanism of the migration- 
proliferation dichotomy was introduced for a comb model 
[1 81 ] , which is a toy model, where the migration prolifera- 
tion dichotomy is naturally apparent. The collective be- 
havior of cells was considered and the primary focus was 
on the influence of cell fission on the transport properties 
of cells. 



II. CELL KINETICS IN PRESENCE OF THE 
TTF 

In the present study, we first consider the kinetic prop- 
erties of cancer development in the presence of the TTF. 
These can be described by fractional kinetics in the 
framework of the comb model, where it is easier to draw 
an intelligible picture of interplay between high-motility 
aggressive cancer cells and the TTF in the outer-invasive 
region. Following toy model consideration [l8J, we re- 
strict our model to one dimensional subdiffusion in the 
x direction. In the comb model, this anomalous diffu- 
sion can be described by the 2D distribution function 
P = P(x,y,t), and a special behavior is that the dis- 
placement in the x-direction is possible only along the 
structure axis (z-axis at y = 0). The Fokker-Planck 



equation in some dimensionless variables reads 

d t P = 6(y)d 2 x P + dd 2 y P , 

where d is an effective diffusion coefficient. Obviously, 
cell division is random in x and discontinuous, contrary 
to that in the tumor core. It can be reasonably considered 
as a fractal set F v {x) with the fractal dimension v, which 
is embedded in the ID space, < v < 1. Therefore, 
the effective diffusion coefficient becomes inhomogeneous 
d —¥ dx(x), where x( x ) is a characteristic function of the 
fractal, such that x( x ) = I for x <E F v {x) and \{x) = 
for x ^ F v (x). Now we take into account the influence 
of the TTF that affects (destroys) only the proliferation 
phenotype according to Refs. |l]-u|- Mathematically, this 
process is expressed by diffusion in the y direction with 
decay: 

dd 2 y P(x, y, t) => X (x)[dd 2 y - C]P(x, y, t) , 

where coefficient C defines difference between the prolif- 
eration and the degradation rate. Taking this into ac- 
count, one arrives at a model that makes it possible to 
understand the influence of the TTF on the cancer de- 
velopment 



d t P = 8{y)d 2 x P + X (x){dd 2 v - C]P(x, y, t) 



(1) 



First, we apply the Fourier transform to Eq. ([T]) with re- 
spect to the x coordinate. To apply this transformation 
to the last term in Eq. ([T]), we use the auxiliary identity 

x(x)f(x) = o x J^ oo x(y)f(y)dy = -d x J™x(y)f(y)dy 

with the boundary conditions P(x — ±oo) = 0. This in- 
tegration with the characteristic function can be carried 
out by means of a convolution (L9l - |2l| [361 ] . To this end we 
will use the terminology and useful notations of fractional 
integration and differentiation [T7|, I22H26J . Fractional in- 
tegration of the order of v is defined by the operator (see 
also Appendix) 



>£/(*) 



I» J- 



f{y){x-yy- l dy., 



(2) 



where < v < 1 and T(v) is the Gamma function. The 
Weyl fractional derivative is 

Wl~ V = cU-oo^/W] = -ooI V Ad x f{x)] . 

Using this notation, one introduces the coarse-graining 
integration with the characteristic function in the form 
of the Weyl fractional derivative [19| 

X (x)P(x,y,t)=>W l x - v P(x,y,t). (3) 

The Fourier transform of the Weyl derivative yields 

P x Wl- v P(x, y, t) = (ikf-'Pik, y, t) , 

where IF x \P(x,y,t)\ = P(k,y,t). Presenting this Fourier 
transform in the symmetrical form, one obtains 



d t P 



-k 2 S(y)P- 



[dd 2 P 



CP] 



(4) 



The last term in the r.h.s. of Eq. ((4]) is eliminated by 
substituting P(k,y,t) = e~ c l fe l f "F(k,y : t). Then one 
carries out the Laplace transform in the time domain 
C[F(k, y, t)] = F(k, y, s). Looking for the solution of the 
Laplace image in the form 



F(k,y,s) = cxp[-\y\^/\k\^s/d}f(k, i 



(5) 



one arrives at the intermediate expression in the form of 
the Laplace and Fourier inversions 



P{x,y,t)=f^le 
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( 6 ) 

One has to recognize that the y axis is the auxiliary co- 
ordinate, which determines the cell proliferating process 
(cell fission). Therefore to find the complete distribution 
of cancer cells in the x axis, integration over y is per- 
formed: P(x,t) = J_ P(x,y,t)dy. Both the integration 
over y and the inverse Laplace transform are carried out 
exactly. This, eventually, yields a solution in the convo- 
lution form 



P(x,t) 

~ 2ir J—o 



(cty 



e lkx e~ 



r n 

J — oo 



•f^r-W (7) 

-c\kf-t Ei _ { _ }.^\k^+n/d s jdk . (8) 



Here 



y ' 2ni L u c 



is the Mittag-Leffier function defined by the inverse 
Laplace transform with a corresponding deformation of 
the contour of the integration [23], while V{z) is a so- 
lution for the untreated glioma with the scaling vari- 
able z = x/t^ and K(z) is the kernel of the TTF 
treatment. For the small argument, which corresponds 
(for a short time) to a long-scale tail of the distribu- 
tion, the Mittag-Leffler function decays exponentially 

exp ( — K v y/\k\ 3+u t\ [Lj, I13] with the generalized trans- 
port coefficient K v = [2r(3/2)\/d]- 1 . This yields the 
solution for the untreated cancer in the form of the Fox 
functions [17| . Its large argument asymptotic solution 
yields the power law decay V{z) ~ \/t/\x\ 5+l/ that is 
an indication of the non-localized diffusive cancer in the 
outer-invasive region[37]. 

Now we turn to the non-invasive treatment in Eqs. ([?]) 
and (jHJ, described by the kernel 
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(cty 
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(cty 



= T: 



-C\k\ x -"t 



which is the Fox function. Therefore, the resulting so- 
lution in Eqs. and © is the convolution of two Fox 
functions, which means that the TTF is a fairly ineffi- 
cient treatment in the outer-invasive region. Our main 



aim here is to understand the efficiency of the TTF. As 
seen from Eqs. an< ^ ©, the crucial point is the fractal 
dimension v < 1. For v — 1 the solution is 

P{x,t) = e- Ct V(x 2 /ti), 

where V(z) is the stretched exponent. This means that 
the TTF is most efficient in the multi-cell tumor core and 
in the close vicinity of the core, where the proliferating 
cell distribution is continuous with the fractal dimension 
v = 1. Therefore, to increase the efficiency of the TTF 
in the outer-invasive region, the TTF must act on both 
the migrating and the proliferating cells. 



III. ELECTROPORATION 

In what follows we suggest a possible mechanism of 
such a treatment of the diffusive cancer using the TTF 
in the outer-invasive region. The idea is based on the 
differences between the dielectric properties of normal 
tissue cells and cancer cells that can lead to the local en- 
hancement of the electric field and to electroporation of 
the cancer cells. Therefore, we consider the frequency- 
dependent permittivities of migrating cancer cells e m and 
the normal tissue cells e n . Under certain frequency, the 
condition e rn < e n can be fulfilled. This was found for 
dielectric properties of normal, and malignant lympho- 
cytes [14| . Therefore, the outer-invasive region is a frac- 
tal composite of two dielectrics with permittivities e n and 
e m , where e n = s n (uj) corresponds to the host dielectric, 
while e m = e m (uj) corresponds to the fractal cancer in- 
clusion. We do not consider proliferating cancer cells, 
since the FFT destroys them effectively. Note also the 
wavelength of the TTF is much larger than any brain in- 
homogeneities, and therefore the quasi-electrostatic con- 
sideration is valid. The host volume of normal tissue cells 
is 3D, while the volume of cancer cells is fractal V ~ r fr 
with fractal dimension Df T < 3 due to their low concen- 
tration, which is much less than their concentration in a 
solid tumor. Therefore, the outer-invasive region can be 
considered as a random fractal cancer with an averaged 
fractal density p(r) ~ r * -3 , as e.g., for low-grade as- 
trocytomas. Note, in electrostatics, p(r) is considered as 
a fractal quenched disorder distribution. 

Electrostatics in the frequency domain (when E = 
E(w,r)), in the random fractal dielectric composite is 
described by the Maxwell equation 



Ee(r) 



0. 



(9) 



where the inhomogeneous permittivity is a function of 
the frequency of external TTF and coordinates: e(r) = 
e(r, u). The boundary conditions are determined by the 
external TTF. We work with dimensionless coordinates, 
which are scaled by the cell size. It is convenient to split 
the electric field inside the brain into two components 
E(r) = E(r) + Ei(r), where E is the homogeneous elec- 
tric field with the condition V • E = 0, while Ei is the 



electric field due to fractal cancer inhomogeneities. We 
take into account that the fractal mass M(r) ~ r fr 
results from averaging over all the possible realizations 
of the random fractal. The averaged density distribu- 
tion p(r) is isotropic and depends on the radius only (in 
complete agreement with the probabilistic sense of ran- 
dom fractals [28|, [29| )■ This also supposes that the di- 
electric properties of the random fractal composite are 
isotropic, such that the space-dependent averaged po- 
larization is a function of the radius only: e(r) = s(r). 
Therefore, it is reasonable to suppose that the response 
to the TTF due to the fractal inhomogeneities changes 
as a function of the radius Ei(w,r) = Ei(w,r), as well. 
Thus the divergence in the Maxwell equation for Ei 
takes into account only the r component of the elec- 
tric field. This yields V ■ E x = V r £a(r), where Ei(r) 
is the r component of the electric field (the r index is 
omitted). The inhomogeneous permittivity of the com- 
posite we define by means of the fractal characteristic 
function x( r ) = 1, if T" belongs to cancer cells, and 
x(r) = 0, if r belongs to the normal tissues. Therefore 
e(r) = e m x(r) + [1 - x(r)]£ n = e„[£x(r) + 1 ], where we 
introduced dimensionless parameter £ = (s m — e n )/s n . 
Substituting this in the Maxwell equation (j9)), one ob- 
tains 



[&(r) + l]V r £i + E lX '(r) = -E rX '(r) 



(10) 



where x'( r ) — :r~x( r )> an d we take into account that 



V r = -pi-^fr 2 . We introduce a new function B(r) = 
r 2 E\(r), and then the integration of Eq. (fT0| with the 
elementary volume dV = 4irr 2 dr reads 



£ / x(r)B'dr 
/o 



B + f / B X '{r)dr = 



£ / E r (ry X '(r)dr, (11) 



where B'(r) = ^B{r). 

Now we in a position to use fractional calculus for inte- 
gration with the fractal characteristic function. By com- 
plete analogy with Eqs. @ and ([3]), this integration cor- 
responds to the convolution integral with the averaged 
density p(r) that yields [19J j 



X(r)f{r)dr 



T(D h - 2) 

oI r D "- 2 /(r) 



(r - r') Db - 3 f{r')dr' 



where f(r) is arbitrary, in particular, f(r) — B'(r). The 
situation with polarization terms is more complicated, 
and one ought to take into account the electrostatic prob- 
lem of polarization of dielectric balls by an external elec- 
tric field 0L All details of the inferring can be found in 
Refs. 113,1301, and the finite result reads 



X \r)B{r)dr ^ Po I^'- 2 B{r) . 



where p = jtj is a polarization parameter that is well 
known in the literature 0, |3l|, [32|. Using this expres- 
sion and taking into account that in homogeneous media, 
when cancer is absent, E = Ea/e n , where Eq — Eq{uj) is 
the amplitude of the TTF in the frequency domain, one 
easily calculates the term in the r.h.s. of Eq. (fTTj) : 



pE i 



o/ Av-2 r 2 _ 



2pE i 



„D tt 



s n T(D iT + 1) 

Now Eq. ([lip can be rewritten in the form convenient 
for solving 



Z I? tr ~ 2 !3'+B+PoI? tr ~ 2 B = - 



2 P Eo£ 



e n T(D h . + 1) 



r u ^ . (12) 



We consider a case with negative £ and |£| <C 1, which 
is the most important for the electric responce to the 
TTF. In this case one neglects the poI® tr ~ 2 B term. The 
Laplace transform method yields the solution in the form 
of the two-parameter Mittag-Lefner function [27[ 






r l-P r u oc-P e ur du 
~2riL u°-l/|$| ' 



where, in our case, a — 3 — Df r , j3 = 4. For small |£|, the 
large argument asymptotics of the Mittag-Lefner func- 
tion is exponential [27J : £ a ^(z) ~ ^z~ 3 e z , where z — t^t. 
This, finally, yields the solution for the electric field 



E^r)--^ 



Eo |£| 



o ex P , 



(ii). " = 3-£> fr . (13) 



Therefore, the respond electric field can be large enough 
to break the cell membrane. For example, for |£| ~ 0.2, 
Df r = 2.5, and e n ~ 10 2 , the electric field response 
is 10 4 -r 10 5 V/cm, which exerts the irreversible elec- 
troporation 33] due to the external TTF with ampli- 
tude Eq ~ lV/cm. This can be a mechanism for ab- 
lation of cancer cells, which effectively acts on migra- 
tory cancer cells. An important condition for this re- 
alization is e m (uj) < e n (uj). These permittivities were 
observed in time domain dielectric spectroscopy in ex- 
perimental studies of the static and dynamic dielectric 
properties of normal, transformed, and malignant B- and 
T- lymphocytes [lj]- Such experimental studies of the 
glioma cells can not be overestimated. 



IV. CONCLUSION 

We presented two models of treating glioma by alter- 
nating, radio frequency electric field, which is the tumor- 
treating field (TTF) in the presence of the migration- 
proliferation dichotomy of cancer cells. The first model 
considers the treating tumor development on a comb 
model. In the framework of this toy model of the 
migration-proliferation dichotomy, based on the frac- 
tional cell transport, it was possible to estimate the effec- 
tiveness of the TTF treatment in the outer-invasive re- 
gion of the tumor development. We also show that while 



the TTF is highly effective in the multi-cell tumor core, 
its action is ineffective in the presence of the migration- 
proliferation dichotomy. The key reason is the fractal 
structure of the glioma cancer in the outer-invasive re- 
gion. Therefore, another possible mechanism of the dif- 
fusive cancer treatment by the low-intensity alternating 
electric field has been suggested. The idea is based on 
the difference between the dielectric properties of nor- 
mal tissue cells and cancer cells J14l . [L5J . Since the cell 
permittivities are functions of the TTF frequency, it is 
supposed that there are conditions where the permittiv- 
ity of normal cells is larger than that of cancer cells [lj] . 
Therefore, considering the tumor invasion through nor- 
mal tissue as a fractal dielectric composite in the pres- 
ence of a high-frequency electric field, the electrostatic 
Maxwell equation in a fractal medium was considered. 
An analytical solution was obtained in the framework of 
fractional calculus, and an essential enhancement of the 
electric field was obtained. The result depends essen- 
tially on the cancer fractal dimension Df T and the dif- 
ference between the permittivities |£| < 1. In view of 
the toy model ([T]), this treatment leads to the solution 
P(x,t) ~ e~ ct T'(x,t), where C is electroparation, which 
is an effective treatment in the invasive region. 

A key quantity of the cancer treatment is localization 
of the electroporation field E\(r) inside the cancer. There 
is a straightforward analogy with nanoplasmonics (see 
e.g., [3J, [35|), where the electric field enhancement is 
due to a so-called surface-plasmon resonance for a metal- 
dielectric composite, and localized surface plasmon os- 
cillations are charge density oscillations confined to the 
conducting fractal nanostructure. The essential differ- 
ence and similarity should be admitted: this biological, 
cell enhancement of the electric field is not resonant, but 
geometrical due to the fractal cancer structure [301 . In 
this connection, in vitro experiments can be important 
for further understanding the interplay between the TTF 
and the migration-proliferation dichotomy 



V. APPENDIX 

Fractional integration of the order of a is defined by 
the operator 



JZf(x) 



1 



TO) J a 



f(y){x-y) a -"dy 



where a > 0, x > a and T(z) is the Gamma function. 
The fractional derivative is the inverse operator to a I" 
as a D"f(x) = a Ix a an d alx = aD~ a . Its explicit form 

is 



a^x — 



1 



r(-a) J a 



f(y)(x-v)- 1 - a dy. 



For arbitrary a > this integral diverges, and as a result 
of a regularization procedure, there are two alternative 
definitions of a D~ a . For an integer n defined as n — 



1 < a < n, one obtains the Riemann-Liouville fractional 
derivative of the form 

a D% L f{x) = {<rix n )j n x - a m , 

and fractional derivative in the Caputo form 

a D a c f{x)= a i:- a f^{x). 



When a — — oo, the resulting Weyl derivative is 

One also has - 00 Dy / e x = e x This property is convenient 
for the Fourier transform T [W Q /(a;)] = (ik) a f(k), where 
P[f{x)] = f(k). 
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